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O ■ ABSTRACT 

G^ ■ The fastest ejecta of supernova explosions propagate as a precursor to the main 

supernova shock wave, and can be quite energetic. The spectrum of such fast ejecta is 
Ph| estimated based on recent analytic and numerical supernova models, and found to be a 

O ' power law having a cutoff at an energy of order 10 MeV/nucl, the precise value of which 

depends on the supernova mass and energetics. With cutoffs in this range there can be 
d '. significant flux with energies above the thresholds for 7-ray and Li, Be, and B produc- 

tion. These nuclear interactions should inevitably accompany the passage of prompt 
supernova ejecta through the surrounding medium. The fast particle composition is 
that of the outermost layers of the progenitor; if the progenitor experienced significant 
loss of its envelope to a companion or by mass loss, then the composition is nonso- 
lar, and in particular, metal-rich. Such a composition and spectrum of fast particles 
is required to explain the recent COMPTEL observations of 7-ray emission in Orion. 
Supernova ejecta from one progenitor having lost a large amount of mass are shown to 
quantitatively account for the Orion 7 rays, and to imply the presence of other, weaker 
lines and of a supernova explosion in Orion in the last < 10^ years. Implications for 
this mechanism as a mode of nucleosynthesis, and as a potential supernova diagnostic, 
are discussed. Model dependences and uncertainties are noted, and the need is shown 
for accurate measurements of 7-ray and spallation production cross sections at energies 
near threshold. 



Subject headings: acceleration of particles — gamma rays:theory — nuclear reactions, nu- 
cleosynthesis , abundances — sup ernovae : general 
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1. Introduction 

An exploding supernova is one of the most 
powerful engines in the Galaxy. Modern super- 
nova models predict that the ejecta should have 
a distribution in velocity, with the outermost lay- 
ers being ejected at great speed. Indeed, it was 
originally hoped (e.g., Colgate & Johnson 1960 ) 
that these energetic particles might be the cos- 
mic rays, and thus supernovae would be the di- 
rect agents of cosmic ray acceleration. This has 
been shown to be incorrect (e.g., Casse & Goret 
1978| ; Meyer |l985a| , |l985bD , and instead it is now 
believed that the supernova blast leads indirectly 
to cosmic ray production via shock acceleration 
(as reviewed, e.g., in Blanford & Eichler 1987 , 
and Jones & Ellison |T99l] ). 



Nevertheless, very fast ejecta are indeed cre- 
ated by supernova explosions. Indeed, some 
have energies above thresholds for nuclear in- 
teractions, in particular for 7-ray and Li, Be, 
& B (LiBeB) production. These interactions 
should inevitably accompany the passage of fast 
ejecta through the medium surrounding the su- 
pernova. Given the presence of these interac- 
tions, the crucial question is whether they are 
numerous enough to be detectable. Indeed, if 
the flux of the fast particles is sufficiently high, 
the 7-rays could offer a new and unique window 
on the explosion, and the LiBeB could provide a 
record of such events, integrated over the history 
of the Galaxy. 

In this paper we estimate the flux of the fast 
supernova ejecta and find that it is in fact suffi- 
ciently energetic and intense to lead to significant 
7-ray production, and less significant LiBeB pro- 
duction. It is thus plausible that this mechanism 
is of great interest as a diagnostic of supernovae 
in general. More particularly, it is interesting in 
light of recent Orion 7-ray observations, which 
can be accounted for by prompt ejecta a single 
type Ic supernova (if it indeed arises from a CO 
core progenitor as recently suggested by Nomoto 



et al. 1994). 



2. Composition and spectrum of ener- 
getic massive star ejecta 

Consider a supernova exploding in some medium, 
e.g. the ISM or, in the case of Orion, a molecular 
cloud. The ejecta will interact with this medium, 
which may have been altered locally by the su- 
pernova progenitor. For example, if there was a 
strong wind, it will have cleared out a cavity in 
the local medium in which the local material is 
replaced by that of the wind. Moreover, a typi- 
cal progenitor begins life as an O star whose UV 
radiation partially clears out the region near the 
star. In fact, it is only important for our purposes 
that the explosion be in some medium, and we 
will see that the particles we are interested in will 
travel much further than the lengthscale of local 
inhomogeneities . 

After the explosion the bulk of the material 
will form a shock and pass through the usual 
phases (free-coasting, Sedov, etc.). But the pre- 
cursors to the blast, the fastest particles ejected, 
are dilute and will not participate in the shock. 
They will instead propagate out and begin to lose 
energy predominantly through ionization losses, 
eventually to be stopped. These are the particles 
in which we are interested. Specifically, we want 
to calculate the intensity, spectrum, composition, 
and duration of the precursor nuclei. 

Before making detailed calculations, let us fix 
the order of magnitude of the effect. Typical su- 
pernova velocity spectra have maximum veloci- 
ties of the order of 40,000 km/s, i.e., with particle 
energies (per nucleon) up to 8 MeV/nucl. Typi- 
cal 7-ray thresholds (for reactions between an a 
particle and a heavy nucleus) are ~ 2 MeV/nucl, 
and spallation thresholds a few times higher; pro- 
ton thresholds are higher still. A threshold en- 
ergy of E'th = 2 MeV/nucl corresponds to a 
velocity of Vth = 20 000 km/s; the amount of 
mass ejected above this energy is roughly M{v > 
Vth) ~ O.OOIM0, i.e., a tenth of a percent of the 
total ejecta. This material has a kinetic energy of 
order 4 x 10^^ erg, i.e., about 0.4% of the super- 
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nova energy. This is to be compared with a typi- 
cal efficiency of ~ 1% for acceleration of galactic 
cosmic ray supernova shocks (e.g., Gaisser 1990). 



Since prompt ejecta have comparable energetics, 
it is plausible that their effects could be impor- 
tant. 

We are thus interested in calculating the ini- 
tial spectrum of particles and determining their 
nuclear interactions as they propagate away from 
the explosion. The initial spectrum is calculated 
as follows. Imagine as well that the star emits at 
total number Ni of isotopes i. The spectrum is 
obtained from the calculations of the mass, veloc- 
ity, and composition of supernova shells. Specif- 
ically, we know the velocity v{M) and mass frac- 
tion Xi{M) as a function of mass coordinate M. 
We want the number spectrum dNi/dE for each 
species i of interest. We have, for the total num- 
ber of particles. 
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integrated over the ejected shells mass between 
the inner Mi and outer mass coordinates of 
the ejecta. But by definition we also have 
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and so by equations ([l|) and (|2|) we have 
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and in the nonrelativistic case which holds here, 
we may simplify further to 
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of LiBeB or 7-ray production by this flux. For- 
mally, the rate is given by 



dt 



rex. 



(5) 



Now we may define and then estimate the rate 



Note that the spectrum dNi/dE will be out of 
equilibrium and so time dependent (in contrast 
to the steady state of galactic cosmic rays) . Be- 
cause of the energy losses, the fast particle fiux 
will be continuously degraded in intensity and 
energy over time. For the purposes of a first ap- 
proximation, however, imagine a constant flux 
lasting for a duration Ai. The appropriateness 
of this approximation will be discussed in §^ 

From eq. (P) it is clear that what is crucial 
to our calculation is the nature of the fastest 
ejecta, i.e., the composition Xi and the veloc- 
ity spectrum dM / dv of the outermost supernova 
shells. Both of these will be sensitive to the pro- 
genitor mass and history. Speciflcally, the outer- 
most composition (and indeed the pre-explosion 
mass) will depend on the amount of mass loss 
the progenitor has suffered. Although the out- 
ermost layers will be born with the composition 
of the ISM at the time, with sufficient mass loss 
this could be shed and the inner nucleosynthetic 
shells can be exposed. Thus mass loss considera- 
tions are crucial for determining the fast particle 
composition. 

Since mass loss can change the progenitor 
mass considerably, it also affects the velocity of 
the ejecta, and thus influences the normalization 
of the velocity spectrum. However, it appears 
that the shape of the velocity spectrum is gener- 
ically a power law to a reasonable approxima- 
tion. In a beautiful analytic treatment, Nady- 
ozhin (1994) has examined the passage of the 
type II supernova shock through the outermost 
layers. He shows that in this a self-similar solu- 
tion holds, and calculates the velocity spectrum 
of the fastest layers. He finds the spectrum to 
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have the form of a power law with a cutoff at a 
maximum velocity fmax- Specifically, he finds 



M{v) 



j B{V " - Umax + AM , V < i;max 
\ , V> -Umax 

(6) 

where ;B is a constant that depends on the details 
of the hydrodynamics, while Vmax scales roughly 
as Umax ^Ju /Me] with the supernova mechan- 
ical energy U and the ejected mass MejJ^ For 
SN 1987A, Nadyozhin calculates fmax = 

(26 - 

40) X 10^ cm/s, corresponding to an energy (per 
nucleon) cutoff of -Emax = 6 — 8 MeV/nucl. 

As we will see in detail below, the presence 
of the cutoff velocity Umax is crucial in determin- 
ing the nuclear interactions of the fast particles. 
Physically, the cutoff arises from the quenching of 
the shock's acceleration of the outer layers of the 
supernova. As the shock moves through the de- 
creasing density of the outer atmosphere, the ma- 
terial is accelerated progressively faster in reac- 
tion to the large radiation pressure of the shock. 
This process ceases, however, when the shock ap- 
proaches the surface and the optical depth for the 
shock radiation becomes comparable to the opti- 
cal depth to reach the photosphere. At this point 
the radiation escapes from the putative super- 
nova as an ultraviolet flash, and the shock causes 
no further direct acceleration (for further discus- 



sion see Nadyozhin 1994 ). 

A very nice feature of the Nadyozhin spectrum 
is that it takes the exceedingly simple form of a 
power law. Thus one can immediately deduce 
that the effective particle spectrum (cf. eq. ^ is 
also a power law dM/dv oc v~'^^~^^\ or 

^(i^)oci?-(^+l)/^ i?<i?max (7) 

dv 

For SN 1987A, Nadyozhin estimates that s = 
7.2, which gives an energy dependence of E~^'^. 
This result is supported by Ensman & Burrows 



^ Note there is an additional shell ejected at the highest 
velocity, but its mass is so small (AM ~ IO'^Mq) that 
we ignore it. 



(1992) whose numerical model of SN 1987A was 
specifically designed to calculate shock breakout. 

The slope s of the velocity spectrum (eq. P) 
is related to the hydrodynamic properties of the 
self-similar shock solution. In particular, s = 
(n + l)/7, where n is a power law index for the 
density dependence on the distance to the sur- 
face, and A is a similar index for the velocity 
(related to n and the adiabatic index). At any 
rate, the universality of the power law solution is 
very fortuitous, as it allows us to make general 
remarks about all SNe. 

The other major feature of the velocity spec- 
trum which is equally crucial for our purposes is 
the existence of the cutoff in velocity, and thus in 
energy. The cutoff is quite near to the 7-ray and 
spallation thresholds, and so will greatly empha- 
size the differences among the different thresh- 
olds. To be sure, there is the above-mentioned 
scaling with supernova mass and energy, so a 
lighter pre-supernova star could perhaps achieve 
higher energies. Regardlessly, the low energy be- 
havior of the cross sections will be crucial. In 
particular, protons reactions are disallowed be- 
low ^ 20 MeV/nucl for C and O reactions, but 
since a thresholds are lower by a factor 4 in en- 
ergy per nucleon, the components of the ejecta in 
a particles, and heavier nuclei which can interact 
with interstellar He, will be crucial. 

The smallness of the cutoff energy means that 
the numerical calculations are difficult to do cor- 
rectly. Results are very sensitive to the gener- 
ally unmeasured threshold behavior of the cross 
sections. Semianalytic fits to experimental data 
(Read & Viola 1984) do extrapolate down to 
threshold, but clearly there is a need to under- 
stand this low-energy physics exceedingly well, 
and we urge experimentalists to consider these 
measurements. Indeed if accurate cross sec- 
tions were available, one could turn the problem 
around, and use the 7-ray and spallation yields 
as a powerful diagnostic on the fast ejecta's spec- 
trum and composition. 

Another word of caution: Nadyozhin's calcu- 
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lation was plainly not intended for the use we 
make of it. In particular, it is doubtful that the 
real cutoff behavior is as sharp as he calculates 
it to be (eq. ^). In actual supernovae it is likely 
that the quenching of the shock is a continuous 
process and so while the cutoff scale may be cor- 
rect, the details of the spectrum at these highest 
energies may be quite complicated. One should 
thus regard numerical results obtained via eq. @ 
as rough indicators of the yields one expects from 
these fast ejecta. Nevertheless, we forge ahead in 
applying our method to the recent 7-ray obser- 
vations of Orion. 

3. Gamma-ray Observations of Orion 

The 7-ray production in the Orion complex 
is huge, as Bloemen et al. ( [1993| ) recently deter- 
mined with the COMPTEL experiment on the 
Compton Gamma Ray Observatory. The most 
surprising feature of the 7 emission from Orion 
is that the large flux is confined to the 3 — 7 MeV 
range. Bloemen et al. ( |1993| ) measured the flux 
in the 3-7 MeV band to be 

$^(3-7 MeV) = (1.0 ±0.15) X 10"'' cm^^ . 

(8) 

In this range, the spectrum shows features con- 
sistent with de-excitation emission from the 4.44 
MeV excited state of ^^C* and the 6.13 MeV state 
of ^^O*; these features have roughly the same am- 
plitude. 

The features at 4 and 6 MeV appear wide, 
which would have implications for the nature of 
the accelerated particles due to the kinematics of 
the nuclear excitations. Fast C and O nuclei give 
rise to fast ^^C* and ^^O* which have Doppler- 
broadened decay lines, whereas fast protons and 
a nuclei lead to slow ^^q* g^j^j 16q* j^aving nar- 
row lines. However, Bloemen et al. ( p.993| ) cau- 
tion that their lineshapes are not corrected for 
broadening due to finite detector response, and 
to the widths are artificial, a point also noted by 
Reeves & Prantzos ( |1995D and emphasized via 
detailed calculations by Ramaty, Kozlovsky, & 



Lingenfelter ( |l995a| , p95b| ). 

The COMPTEL experiment also set limits to 
the flux outside of the 3 — 7 MeV band; in par- 
ticular, the flux in the 1 — 3 MeV range was con- 
strained to be < 1/7 that of the 3 — 7 MeV range. 
The lack of such lines is a crucial indicator of the 
lack of fast protons or a particles, which would 
excite ambient Fe, Si, and Ne, leading to emission 
in this range. 

It is assumed that these lines have their ori- 
gin in energetic particle interactions with ambi- 
ent matter in the Orion complex, but it is clear 
that any explanation of the 7 emission will have 
to be unconventional. No lines are observed be- 
low the 3 — 7 MeV band, and no excess vr'' decays 
have been observed at higher energies (i.e., no ex- 
cess over the level expected from galactic cosmic 
rays; see Digel, Hunter, & Mukherjee |l99^ ). The 
if galactic cosmic rays produced the 3 — 7 MeV 
lines, they would also produce 1-3 MeV 7-ray 
lines which are not observed; also the GCR inten- 
sity would have to be much larger, which would 
give a high excess of vr'' decay photons, again un- 
observed. Moreover, the energetic demand would 
be enormous. This process is therefore excluded. 

The most economic way to obtain the ob- 
served 7-ray flux is through accelerated C and 
O nuclei, and above all, a suppression of pro- 



tons (Ramaty, Kozlovsky, & Lingenfelter 1995a , 
1995b| ; Vangioni-F lam e t al. |l99|; C asse, Lehoucq, 



& Vangioni-Flam |199E| ; and Casse et al. |1994|) 
Indeed, if one assumes merely that the Orion 
7-rays arise from energetic particle interactions, 
then one may deduce in a model-independent 
way the particles must be a large flux with low 
energy, predominantly composed of C and O 
rather than protons and a particles. 

A low energy CO flux was suggested by the 
Bloemen et al. ( |1993| ) paper itself, and the low 
energy character was assumed in all subsequent 
efforts at explaining the 7-ray emission. Bykov 
& Bloemen ( |1994| ) and Nath & Biermann ( |1994|) 



connect the observation to theories of cosmic ray 



origin. Bykov & Bloemen (1994) propose a spe- 
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cific spectrum arising from collisions between su- 
pernova shocks and stellar winds. They use pure 
supernova compositions in their calculations, but 
it remains unclear how much of an admixture of 
normal material will be present in the hot re- 
gions they consider. Nath & Biermann ( |1994| ) 
favor acceleration in shocks of winds of hot O 
and B stars; they argue that this naturally gives 
the desired composition. Ramaty, Kozlovsky, & 



Lingenfelter (1995a) pointed out that while the 
observed 7-ray features appear broad, there is 
insufficient experimental resolution to exclude a 
narrow line component. Parameterizing a low en- 
ergy spectrum, they tried different source com- 
positions, showing that a composition like that of 
a WC supernova wind would be needed to meet 
the 7 spectrum and energetic constraints. These 
authors have also suggested (Ramaty, Kozlovsky, 
& Lingenfelter 1995b| ) that the low energy cosmic 
rays could have an origin akin to that of the so- 
called anomalous component of the cosmic rays 
observed at earth (a suggestion also made by Jin 
h Clayton |l99^ ). 

A low energy cosmic ray flux has implica- 
tions for the nucleosynthesis of Li, Be, and B 
(LiBeB), as first pointed out by Meneguzzi & 
Reeve s (|1975b| ). Casse, Lehoucq, & Vangioni- 
Flam (1995) investigated the implications of such 
a process for LiBeB chemical evolution. They 
show that the LiBeB synthesis can be significant 
(for an exposure time of 10^ yr), comparable to 
that of galactic cosmic rays (e.g. Walker, Math- 



ews, & Viola [19851) . 

Much of the work on the Orion 7-rays thus 
far has either: (1) assumed arbitrary flux shape, 
or (2) assumed SN composition. These first ap- 
proaches adopt a spectrum of the form oc , 
with s = below an energy E^, and n = 5 — 10 
above (a form chosen for exploratory pur- 



poses, following Meneguzzi & Reeves (1975b) and 
Ramaty, Kozlovsky, & Lingenfelter ( 1979| )) and 
try different supernova compositions, as Bykov 
& Bloemen (1994). We propose a more specific 
scenario related to a type Ic supernova explosion 



in a molecular cloud. 

Our mechanism naturally has all of the needed 
features, provided only that the supernova ex- 
ploding within Orion be sufficiently massive to 
have undergone significant mass loss. Then the 
progenitor would be stripped to a CO core, which 
would be ejected at high speeds (for a fixed ex- 
plosion energy, a smaller core mass has less grav- 
itational binding energy and so leads to higher 
kinetic energy for the ejecta). Furthermore, nei- 
ther the fast particle spectrum nor the number of 
particles is arbitrary, but rather these are given 
by the supernova model. Thus the predictions of 
this process are specific and so the mechanism is 
testable. 

4. Confrontation with the Orion Obser- 
vation 

We have seen that theory allows us to compute 
N^, while the Orion observation report the flux 
<I>-y at Earth. These quantities are simply related 
by 

dK 



= At: 
dt ^ 



39 „-l 



(3.0 ±0.45) X 10^^ s 



(9) 



using i?orion ~ 500 pc as the distance to Orion, 
and assuming the flux to be emitted isotropi- 
cally so that the total subtended solid angle is 
47r. Of course it is conceivable that the emission 
is anisotropic due to, e.g., magnetic flelds hav- 
ing somehow focussed (or defocussed) the beam 
in our direction. In the lack of clear evidence for 
such effects, we will henceforth make the simplest 
assumption regarding the angular distribution of 
the radiation and consider it to be isotropic. 

To get a feel for the numerical results, let us 
make an order of magnitude estimate of the 7- 
ray emissivity. The particles of interest are the 
ejecta having a velocity larger than the reaction 
thresholds; these particles have mass Mej{v > 
Vth)- Taking an energy (-Emin = 2 MeV/nucl) 
slightly above the threshold for gamma rays 
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1.5 MeV for a+C), we have 



N, 



7 ' - norion V Cr'^c ^CO {v > Vmin) 

^ X^O Meiiv>Vmin) 
= J^Orion 0"pC - 



1.0 X 10^^ Xc 

'^Orion 



(10) 



n 



20 cm-^ 10 mb 0.03 Mq 

where we have taken the cross section to be a 
typical one in the low energy range, about 10% 
of the maximum, and have allowed for an aver- 
age density in Orion that is moderately enhanced 
over that of the ISM. 

Note that eq. ( |lO|) shows that for the fiducial 
values we have chosen, the emissivity is just at 
the observed levels for the C lines (which con- 
tribute about half of the total emissivity) but 
only if Xq (and similarly Xq) is close to unity. 
That this, the fast ejecta must be predominantly 
composed of C and O. But the fastest ejecta 
arises from material at the outermost layers of 
the supernova progenitor. One expects a "nor- 
mal" supernova progenitor, one not having sig- 
nificant mass loss, to have a composition similar 
to that of the local ISM at the time of the pro- 
genitor's birth. This would give C and O mass 
fractions < 10~^, which leads, via eq. (|To[), to 
a significant underproduction of ^^q* g^j^^j 16q* 
lines. 

However, if the progenitor has suffered signifi- 
cant mass loss, or has lost its outer envelope to a 
companion star, then it is possible that the inner 
nucleosynthetic shells could be left at the outer 
layers of the progenitor. In this case the fastest 
ejecta would have a very nonsolar composition. 
In particular, if the progenitor is reduced to a 
CO core, then it would have Xc ^ Xq ~ 0.5; 
this is just the composition required to explain 
the Orion observations. While these nuclei dom- 
inate, there remain some a particles and also a 
small amout of Ne (having the form ^^Ne rather 
than 20Ne). 

We thus suggest that the Orion 7-rays are the 



result of an explosion of a CO core supernova, 
similar to that observed recently in the type 
Ic supernova 19941 (Schmidt, Challis, &: Kirch- 
Clocchiatti, Brotherton, Harkness, & 



1994 



19941 ). Nomoto et al. ( |1994| ) argued that 



ner 

Wheeler 

the observed characteristics of this event are best 
understood in the context of a CO model arising 
from envelope loss to a companion (as opposed 
to mass loss of a single star). While such events 
are not the most common fates of massive stars, 
they are not so unlikely: Nomoto et al. estimate 
a ratio of type Ic to type II explosions of ~ 10%. 

We now apply our formalism in more detail 
to the particular case of a CO core explosion. 
We use the results from the supernova model 
of Nomoto et al. ( p94D (specifically the v{M) 
and Xi(M) relations) to generate the spectrum 
shown in Figure |l|. As expected, most of the par- 
ticles indeed lie at low energies, and the compo- 
sition varies over energies reflecting the variation 
of the composition with mass shell. Since the 
reaction thresholds are all > 2 MeV/nucl, the 
only part of the flux relevant to our problem is 
that at and above this scale. At these energies, 
the composition does not vary with radius (being 
that of the outmost shells). Moreover, the spec- 
trum has, fortunately, settled down to a fairly 
smooth shape close to a power law, with a slope 
of approximately ~ 4, in reasonable agreement 
with the Nadyozhin analytical calculation. Pre- 
liminary calculations using the velocity profiles 
for Wolf-Rayet supernovae (Woosley, Langer, & 
Weaver |l99^ ) give a similar slope and so give fur- 
ther credence to the analytic solutions. 

To investigate the effect of the flux, it is cru- 
cial to know its behavior for energies above those 
specified in the Nomoto calculation. We have 
chosen to extrapolate the flux to higher energies 
by taking a power law with the slope of the last 
points. Specifically, if we write (j) ~ E~°', we 
find a ^ 3.5 around 1 MeV, steepening to a ~ 4 
around 2 MeV; it is this last slope we used in the 
extrapolation. 

Following Nadyozhin, we wish to introduce a 
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maximum energy -Bmax to the flux. The self- 
similar solution provides the scaling mentioned 
above, -Emax oc (f//Moj)-^/^, which Nadyozhin 
normalizes to SN1987A hydrodynamic models. 
Scaling his result to our CO core explosion, and 
assuming a total energy U = 10^^ erg, gives 
Eraax ^ 40 MeV/nucl. However, since the details 
of determining .Bmax are difficult, the supernova 
energy is uncertain as well, we will take -Emaxto 
be a parameter and we will investigate values in 
the range 8 MeV/nucl < E^^ax < 50 MeV/nucl. 

Using this extrapolated flux we have evaluated 
numerically the ^^c*, ^<^0*, 22Ne*, and LiBeB 
production. We assume that the escape length 
Aesc is large and so the LiBeB produced are all 
either thermalized through ionization losses (the 
dominant process), or lost through inelastic col- 
lisions. In fact, since the ionization losses dom- 
inate, the results are insensitive to the escape 
length. Further we assume that the medium has 
a solar system composition as given by Anders & 
Grevesse ( |T989| ). 

Given these parameters, the ratios among the 
7 and spallation rates are fixed; however, to com- 
pute absolute rates requires specification of the 
Orion H density nn- Specifically, we find, for the 
total 7 production from ^^q* _|_ 16q*^ have 

^ = 2.4 X lO^s photons s-^ (11) 

dt 1 cm-3 ^ ^ ^ 

with very little variation over the chosen range 
of -Emax- Comparing this to the measured value 
(eq. P), we see that our mechanism can account 
for the observed Orion 7-ray flux if the ejecta 
encounter an average density 

nn ^ 12 cm~3 . (12) 

This is to be compared with to the estimated av- 
erage density of Orion, n^™'^ ~ 5 cm~^ (Goudis 
1982| ), the uncertainty in which is large enough 
so that it is in good agreement with eq. (|l2|). We 
note, however, that the density requirement rises 
with the time since the explosion, since the flux is 
constantly decreasing in intensity. As discussed 



below (§^) the timescale for significant loss is 
energy-dependent but at least of order 1 kyr. 

Figure §(a) shows calculated the emissivity 
from excited ^^C, ^^O, as well as ^^Ne, plot- 
ted function of ii^max- The normalization 
is fixed by requiring the total ^^q* _|_ 16q* gj^^is- 
sivity to match the Orion observations (eq. ^); as 
discussed in the previous paragraph, this fixes a 
value of riH for each -Emaxj which appears in Fig- 
ure |2|(b). We see that over the preferred range of 
-Emax, the 7-line ratios, and the required density, 
show almost no variation, and what little change 
there is has settled down by 10 MeV/nucl. We 
find C*/0* = 2.7, consistent with the Bloemen 
et al. (|1993D observations. 

It is insufficient to show that the flux can pro- 
duce the observed 3 — 7 MeV lines; we must also 
show there will not be lines outside this range. 
In particular, we must show that the 1 — 3 MeV 
lines are below the COMPTEL hmits (c.f. §|]). 
Here again, the composition of the stripped su- 
pernova core proves advantageous. In particular, 
not only are the fast ejecta dominant in C and O, 
as discussed above, but they also have relatively 
small abundances of heavier elements. In partic- 
ular, the next most abundant element, after C, O, 
and He, is ^^Ne. Although its abundance is not 
large, it is much larger than its solar proportion, 
indeed much larger than the solar proportion of 
heavy elements. Thus this will be the most im- 
portant contributor of heavy element lines, even 
more than those caused by the a reactions in the 
ISM. 

Indeed, we find (C* + 0*)/Ne* > 13 for all 
values of i^max, which falls safely below the ob- 
servational limits, which are of order 7. However, 
if ours is the mechanism for making the C and O 
7-rays, then there should be ^^Ne lines present 
at a level just below the current limits. While 
COMPTEL may not have a good chance to see 
them, INTEGRAL will see them in detail. This 
will then provide a clean discrimination between 
our idea and others, as our prediction of ^^Ne 
lines (at 1.3 MeV), rather than ■^''Ne lines (at 1.6 
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MeV), is likely to be unique. 

Given the roughness of our estimate, and its 
closeness to the observational limit, it is clear 
that the 1 — 3 MeV window offers an important 
constraint. More detailed modeling is needed and 
is under way (Ramaty, Vangioni-Flam, Fields, & 
Casse 1996| ). Indeed, we see that the generic sen- 
sitivity of the 7-ray emission to the supernova 
composition and velocity structure makes the 7 
spectrum a powerful diagnostic tool. Since all 
supernovae should have fast ejecta, one could in 
general turn the problem around and use the con- 
straint as a probe of supernova ejecta. 

Although a full study of fast particle inter- 
actions from all possible supernovae is beyond 
the scope of this paper, nevertheless some dis- 
cussion of model dependence is in order. We are 
interested in the sensitivity to progenitor mass 
and mass loss, as manifested in the flux inten- 
sity, spectrum, and composition. To begin to 
address this question we have studied variants to 
the supernova model we have used thus far. 

The hydrodynamics of the explosion deter- 
mines the intensity and spectrum, and the Nady- 
ozhin ( [1994| ) model suggests (eq. ^ that the 
spectral slope is fairly independent of the ini- 
tial mass, but the velocity scaling and energy 
cutoff depends on Wmax oc ^Ju /Mey The scal- 
ing with mass means that significant mass loss 
can be important since, for a fixed explosion en- 
ergy, it reduces the gravitational potential to be 
overcome and so increases the final kinetic en- 
ergy. To test the accuracy of the expected scal- 
ing with the explosion energy (at a fixed mass), 
we compare the spectrum for a 2.1Mq core at 
energies of (0.6,0.8, 1) X 10^^ erg (Nomoto et al. 
1995a , 1995bD . Results appear in figure | 



with progenitor mass (for a fixed explosion en- 
ergy), we compare the spectrum from the 2.1Mq 
CO core we have considered so far with that of 



we see 



that while the scaling is not exact, it is a good 
approximation. Note that changes in energy scal- 
ing do not affect the intensity of the scaled fiux, 
but the shifts in the velocity scale change the in- 
tensity of the relevant, highest energy fiux. 
To test the accuracy of the expected scaling 



CO cores with masses 1.8 (Iwamoto et al. 1994) 
and 2.9Mq (Nomoto et al. |1995a|) . The results 
appear in figure ^ Here we find the agreement 
is not as good; indeed it appears that the parti- 
cle fluxes are more similar without the scaling]^ 
While for our calculations we have used the nu- 
merical results anyway, it would be interesting 
if the results were less sensitive to ejecta mass 
than expected, as this would lead to larger fluxes 
(above threshold) than expected. At any rate 
we find that clS 3f whole, the slope, the energy 
cutoff, and the intensity is well-described by the 
very simple Nadyozhin model itself coming from 
well- understood physics. 

We now turn to the composition, which de- 
pends strongly on the mass loss, and so is tied 
to much less certain physics. Indeed, mass loss is 
poorly understood and simply treated paramet- 
rically. We have tried using different mass loss 
rates; results appear in Table At any rate, it 
is apparent from figure |^, the composition of the 
fastest shells is fairly homogeneous. Thus one 
need only determine the (uncertain) outermost 
composition without having to get the full pro- 
file and convolving it with the hydrodynamics. 

5. Propagation and Energetics 

Recall that in calculating the rates for nuclear 
interactions of the fast particles we have assumed 
the fiux to be present at full strength for some 
time interval At, during which the spectrum is 
constant. While this is adequate for estimates of 
emissivity, it is too simplistic. At the energies of 
interest the particles are stopped predominantly 
by ionization energy losses. As these losses are a 
strong function of energy and of A and Z (scaling 



^Note, however, that the plot only demonstrates the scaling 
of the overall intensity; the models on which the calcula- 
tions are based do not include the physics of the shock 
breakout and so cannot show the energy cutoff. 
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as Z'^/A), the particle spectrum and composition 
will change with time. 

The calculation is done in detail in Fields 
(|1996| ). The basic effect is that the spectrum con- 
tinually becomes less intense, but also harder. To 
see this, imagine an (unphysical) spectrum that 
is a (5-function in energy. It would degrade from 
its initial energy according to the usual Bethe- 
Bloch dE/dt formula. Thus for a continuous 
spectrum, one expects the entire flux to degrade 
in energy. But since dE/dt oc E~^ , lower en- 
ergy particles lose energy faster, so the spectrum, 
though less intense, has an increasing average 
energy per particle. 

An important related point is that the stop- 
ping ranges are larger than size of Orion. This 
means that if the fast particles are still in Orion 
then the supernova is young. More quantita- 
tively, the range of a 2 MeV/nucl oxygen nucleus 
in Orion is about 50 pc; it traverses this distance 
in about 2 kyr; this is about a factor 3 smaller 
than the crossing time, but it is also a lower 
bound; for particles above about 5 MeV/nucl 
the range is larger than Orion. Thus we can use 
the crossing time to put an approximate upper 
limit to the supernova date of ~ 10^ yr ago. Of 
course, an accurate calculation must include the 
fact that the exposure time is not the total stop- 
ping time but rather the time it takes the fast 
particles to pass below the threshold for 7-ray 
production. Since the stopping time is a strong 
function of energy (t oc £^^) this difference will 
only be important for the particles very close to 
the threshold energy. This effect is discussed in 
Fields (|199|). 

The possibility of a relatively recent supernova 
is constrained by observations of short-lived ra- 
dioactivities. In particular, the lack of observed 
^^Ti (r = 78.2 yr) sets a lower limit to the 
supernova age of order r, assuming that there 
is significant "^^Ti production in the supernova 
in question. In fact, we can make this point 
more quantitative. The CO core supernova mod- 
els we have been using give ^"^Ti yields of order 



M44 = 5 X 1O~^M0 - 2 X lO-^M© (Nomoto et 



al. |1991| ; Kumagai et al. |1991| ; Thielemann et al 



1996), a value close to the 10 "^Mq inferred from 



observations in Gas A (lyudin et al. 1994). If the 
time since the explosion is At, there will be a 
"^^Ti emissivity of 

1 M44 



44 



-At/r 



T 44 TTT-r 



(13) 



at an energy E^ 



1.16 MeV. The observa- 



tional limit on such lines in Orion is L < 3 x 
10^^ photons s^^ (namely the limit for the 1-3 
MeV band), which gives 

M44 



At > rln. 



570 yr 



(14) 



This is a strong constraint, being only a factor 
of 4 shorter than the minimum particle stopping 
time. Taking this result at face value, we are led 
to infer a possible age range of 0.6-10 kyr for the 
supernova. 

However, caution is merited in using the ^*^Ti 
lower bound on the age. The ^^Ti yield is very 
uncertain, as it requires an accurate knowledge 
of which of the deepest layers of the star are 
ejected, and which fall back into the neutron star 
(i.e., the mass cut). Different mass cut choices 
are allowed and can lead to smaller ^^Ti yields 
(although the age limit depends only logarith- 
mically on the mass yield). Also, we note that 
the observational limit we have used is that of 
COMPTEL and is more suited for broad lines; 
Bloemen et al. ( 1993| ) caution that narrow lines 
may be left undetected. And in any case, if the 
age is closer to the ~ 10^ yr crossing time, then 
there is no ^^Ti to be seen by either COMP- 
TEL or INTEGRAL. On the other hand, to take 
more optimistic viewpoint and turn the problem 
around: if there has been a very recent supernova 
in Orion, the ^^Ti lines (or lack thereof) can give 
important information about the mass cut. 

Finally, we turn to the question of energy re- 
quirements. As pointed out by Ramaty, Ko- 
zlovsky, & Lingenfelter ( |19954 |1995b| ), the sig- 
nificant low-energy flux in Orion is a challenge 
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to explain in a way that is energetically econom- 
ical. These authors note that a flux lasting long 
enough (> 10^ yr) to produce significant amounts 
of Li, Be, and B must be composed of relatively 
energetic particles (> 30 MeV/nucl) to avoid an 
unreasonable energy demand. 

Our mechanism also must be shown to have 
reasonable energy requirements. To see most 
simply that indeed the energetics are accept- 
able, recall that our particles have a very defi- 
nite source: one C+0 core supernova. Thus in 
this model the total energy required to acceler- 
ate the particles and thus create the observed 
7-rays is exactly the (mechanical) energy of the 
supernova: 10^^ erg. We note again that the 
flux intensity and composition of the particles is 
fully determined by the supernova model, gener- 
ally and in the speciflc scenario we have adopted 
for Orion. The energy is always fixed to be that 
of the supernova, and so as a consequence we 
are not free to adjust the number of fast parti- 
cles; thus the energy budget is automatically that 
of one supernova, but the resulting fiux may or 
may not be observable (depending notably on the 
ejected mass). 

The number of 7-rays per unit energy injected 
in particles has been calculated by other groups 
and so provides a useful quantity for compari- 
son of the results. In our case, the total energy 
injected in fast particles is above 2 MeV/nucl 
is W = 1.4 X 10^° erg, about 14% of the to- 
tal energy (this is larger than the generic esti- 
mate given in the beginning of §^ due to the 
unusually low ejecta mass). The total number 
of 7-rays made in Orion is the product of the 
present total emissivity and the irradiation time: 
A'^ = N^T ~ lO^'' photons (using a conserva- 
tive exposure time of r = 1 kyr). This gives 
an efficiency N^/W ~ 0.7 photon/erg. By way 
of comparision, Ramaty, Kozlovsky, &; Lingen- 
felter ( 1995a , 1995b| ) find values around 10 pho- 
tons/erg, using a spectrum with a characteristic 
energy of 30 MeV/nucl, which is in fact more 
efficient. 



The implications in this difference in efficien- 
cies can be understood by noting that the fast 
particle propagation is at all not the same in the 
two models. We do not use a steady state flux 
given by a thick target solution, but rater an im- 
pulsive flux. As a result, in our case the expo- 
sure time is fixed as we have just discussed, and 
is in fact much smaller that that of Ramaty et al. 



(I19954 |l995bD or Casse, Lehoucq, & Vangioni- 



Flam (1995). As the flux has a shorter dura- 
tion its total energy requirements are less; on the 
other hand, it also produces less LiBeB, as we 



now see. 



6. Li, Be, and B Production 

With the above discussion one can compute 
the time-integrated production rates. The for- 



malism and discussion is found in Fields ( 1996 ). 
The most important point is that the nuclear in- 
teraction rate proportional to nisM (cf. eq. §), 
but the exposure is set by the stopping time 
which is proportional to tt-ism"^- Thus the pro- 
duction at a given value of energy is related only 
to the stopping power; specifically, it is propor- 
tional to the ionization range as expressed in g 
cm~-^. Roughly speaking, one integrates over ion- 
ization range rather than over time, to arrive at 
a yield that is a given number or mass of LiBeB 
per supernova. 

Thus the absolute yields of LiBeB are indepen- 
dent of the density of the circumstellar medium; 
they are only dependent on the composition of 
the ISM. Consequently, for a fixed composition, 
the yield for a given supernova is independent of 
the region in which it explodes. However, while 
the same argument holds for the total number 
of 7-rays, their intensity does depend on n, as 
we have discussed above (cf. eq. |l^. Note also 
that the composition of the medium, and of the 
ejecta, is crucial. 

We have calculated the LiBeB yields for the 
ISM with a solar system composition. In figure 
we plot the total mass of LiBeB isotopes pro- 
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duced function of -Emax- We see that, unlike 
the results for 7-ray emission, the LiBeB produc- 
tion is very sensitive to i^max- This may be un- 
derstood very simply in terms of cross section 
thresholds. Typically, for spallation processes 
one has ^th(a + N) < Eti,{a + C) < Eti,{a + O), 
so the low -Emax behavior arises from the q + N 
reactions, and the rapid variations happen when 
Ema.x rises above the a thresholds. 

The implications of this mechanism for chem- 
ical evolution of LiBeB are hard to foresee given 
the yields for just one supernova type (and a 
special one at that). Nevertheless, it is of in- 
terest to compare to other suggested LiBeB pro- 
duction mechanisms in supernovae, e.g. the v- 
process yields (Woosley & Weaver 1995| )^. These 
are a function of supernova mass and metallicity, 
but at solar metallicity ^^B production lies in the 
range Mb ~ (1 - 30) x 10"^Mq. It has been 
shown (Olive et al. |1994 ; Timmes, Woosley, & 
Weaver |1995| ) that these yields, when combined 
with the usual GCR production of LiBeB, are 
sufficient to explain the solar ^^B/^'^B ratio (but 
not ^Be). Thus it would seem that since the u- 
process yields are significantly higher than that 
for our particular supernova, it is unlikely that 
our mechanism can produce the solar boron. A 
more detailed analysis of the relative LiBeB con- 
tributions of these processes as well as GCR will 
appear in a forthcoming paper (Vangioni-Flam 
[1996| ). 



et al. 



While the supernova precursors appear to have 
trouble explaining the solar B abundances, it re- 
mains however to be seen if they may explain the 
Pop II abundances. A useful diagnostic is the 
B/0 ratio. For the CO supernova we have con- 
sidered, the O production is Mq ~ 0.2 — OAMq 



In fact, other mechanisms of supernova production of 
LiBeB have been suggested, most notably that of Colgate 
(197J) and Dearborn et al. (1989). Both of these involve 
synthesis within the shock itself, as opposed to our treat- 
ment of the fcist particles' interaction with the ISM. Each 
has also be en sh own to have appar ently fatal problems, 
by Weaver (1976) and Brown et al. (1991), respectively. 



(Thielemann et al. 1994, 1996). This gives, for 
high -Emaxj B/0 ~ 3x 10~^ by number. Compar- 
ing this with the solar ratio (B/O)© = 1.2 x 10~^, 
we see that even if all of the solar O were pro- 
duced by CO supernovae, the most favorable as- 
sumption, we still find B production to fall short 
by a factor of 2-3, again with realistic LiBeB 
production giving smaller numbers. However, in 
Population II, we have B/0 0.1 (B/0)q ~ 
(3 - 6) X lO-'^. Thus we would require only 
~ 5 — 10% of Pop II supernovae to be of this type 
in order to have significant Pop II B production 
by our mechanism; this is to be compared to the 
supernova type Ic/II ratio of ~ 10% (Nomoto et 
al. [199^. 



This suggests that the Pop II B (as well as Li 
and Be) may have arisen from supernova precur- 
sors. If so, then since the fast particles are C and 

0, then there is a natural explanation of the ob- 
served correlation of Be,B oc O in Pop II stars, 

1. e. a linear slope (rather than the quadratic slope 
naively expected from galactic cosmic ray pro- 
duction). This possibility is intriguing and de- 
mands further investigation (Vangioni-Flam et 
al. [1996D. 



Regardless of the impact of our LiBeB yields 
on galactic chemical evolution, the production is 
still high enough to create large local enhance- 
ments of these elements. This is particularly im- 
portant if the fast particles are well-confined by 
local magnetic fields and so the yields are concen- 
trated in a small volume. Indeed, if sufficiently 
localized, the fast particle yields would dominate 
the local LiBeB abundances and ratios. This pos- 
sibility is interesting in light of the Lemoine, Fer- 
let, Sz Vidal-Madjar ( |1993 ) observation of a very 
localized Li isotopic anomaly in the direction of 
( Oph. The observed ^Li/^Li ratio is ~ 2 in one 
of the two observed components; this is tantaliz- 
ingly close to the high E'max values for this ratio 
found in figure ^. If this is so, then the other el- 
emental and isotopic ratios in this region should 
be as given in the figure; in fact, none of these are 
very different from the solar ratios. Thus the sig- 
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nature of supernova precursor production would 
be for the other ratios to appear solar, and of 
course, for the ''Li/^Li measurement to hold up. 

Finally, if the process we discuss does produce 
significant Pop II ^Li, Be, and B, then one con- 
sequence is apparently that there are at least 
two other important sources of these. Galac- 
tic cosmic rays are an inevitable and significant 
source; furthermore, their inability to reproduce 
the ^^B/^'^B ratio suggests the need for another 
mechanism as well. Perhaps the i^-process is the 
remaining source, perhaps some other low-energy 
particle spectrum (Casse, Lehoucq, & Vangioni- 



Flam |1995| ; Ramaty, Kozlovsky, & Lingenfelter 
1995a ), or perhaps something else entirely. At 
any rate, in explaining LiBeB production and 
evolution, one seems to be faced at the moment 
with an embarrassment of riches-more theories 
exist than problems to solve. To resolve the is- 
sue will require not just more observational data 
but also more accurate theory. 

7. Discussion 

We find that the prompt ejecta from a su- 
pernova may be behind the curious 7-ray obser- 
vations of Orion. This mechanism has the ad- 
vantage of providing a definitive physical mech- 
anism for a low-energy particle spectrum with a 
non-solar composition, as seems required by the 
observations. To explain the observation we re- 
quire the progenitor of the responsible supernova 
to have been not a "normal," type II supernova 
with a small wind, but instead an exploding CO 
core star. This leads not only to the observed 
superabundance of low energy CO nuclei, but 
also predicts that there will be ^^Ne (rather than 
^°Ne) lines just below the current upper limits. 
Consequently, our hypothesis is testable. 

Furthermore, via the analytical supernova shock 
models of Nadyozhin ( |1994 ), we see that this 
mechanism is generic to supernova explosions. 
There are always precursor particles preceding 
the bulk of the shock, and these particles in- 



evitably have nuclear interactions with their en- 
vironment if -Emax is sufficiently large. The 7 
production could be a new probe of the super- 
nova explosion, providing a detailed signature 
of the outermost shells of material. This could 
in principle be used as an independent means 
of determining supernova masses (Fields, Casse, 
Nomoto, Schramm, & Vangioni- Flam |199|). 

The implications for LiBeB evolution are un- 
clear at the moment and bear further investi- 
gation. While it is not clear that the LiBeB 
yields will play a major role in galactic chem- 
ical evolution, they nevertheless could account 
for large local variations in these elements. Thus 
this process could offer a potential explanation 
for the low ^Li/^Li ratio observed by Lemoine et 



al. (1993) in the direction of C Oph. 

We caution that calculations we present here 
are rough. Some improvements can be imple- 
mented without too much difficulty. For exam- 
ple, we have treated the propagation of the flux in 
a sketchy manner, but as we indicate the full time 
dependence should be considered. Also, the flux 
first encounters the supernova progenitor's wind, 
and the effect of traversing the wind should be 
included. Both of these effects will be addressed 
in a future paper (Fields |1996 ). 

Another uncertainty in the calculations arises 
from the lack of accurate cross section measure- 
ments for 7-ray and spallation processes occur- 
ring near threshold. We urge experimentalists to 
examine this behavior. Finally, there are some 
uncertainties to our calculations which are more 
difficult to address accurately. The most notable 
of these is the effect of a magnetic field on the 
propagation and stopping of the particles. A siz- 
able field could for example lead to anisotropic 
emission of particles and so of 7-rays. Spiraling 
in magnetic fields will also shorten the total dis- 
placement of the fast particles from the super- 
nova explosion site. These effects are however 
difficult to estimate. 

Despite these problems, the precursor parti- 
cles clearly have the potential to provide a new 
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probe of supernova explosions. More work on 
this mechanism is in order. 
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of the Light Elements). 
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Table 1: Composition of outermost shells, by 
number 



Iniitial Mass 


CO core mass 


Ejected mass 


He 




ieiQ 


^^Ne 


13 Mq 


1.8 Mq 


0.54 Mq 


1.44 


1.39 


= 1 


0.0203 


15 Mq 


2.1 Mq 


0.86 Mq 


1.61 


1.38 


= 1 


0.0337 


18 Mq 


2.9 Mq 


1.54 Mq 


684 


6.42 


= 1 


1.26 
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FIGURE CAPTIONS 



1. The spectrum of the ejecta from a 2.1 Mq 
C+0 core, for a progenitor originally of 
mass 15 Mq. Note the absence of protons 
and the dominance of C and O over a parti- 
cles. The small bounciness is due to numer- 
ical inaccuracy, but the large jumps, e.g. 
in the C spectrum, are real and arise from 
changes in the compositions of the vari- 
ous nuclear burning shells. In the energy 
range of interest {E > 1 MeV/nucl), the 
spectrum is well approximated by a power 
law with a slope of ~ 4, in good agree- 
ment with the analytic results by Nady- 
ozhin ( 1994| ). The endpoints of the curves 
at ~ 3.5MeV are artifacts of the super- 
nova model and do not correspond to i?max- 



2. (a) ^^C*, -^^O*, and ^^Ne* emission as a 
function of -Emax- Curves are normalized 
so that the total CO 7 emission matches 
the Bloemen et al. ( |1993| ) Orion observa- 
tions, i.e., iV(C* + 0*) = 3x 10^9 s'^. This 
amounts to fixing the hydrogen density nn 
at each -E'niax * 

(b) The required hydrogen density nn as a 
function of £'max- Note that the density is 
fairly constant over -Emax! and is only mod- 
estly above the average Orion density. The 
rise at low -Bmax stems from the lower par- 
ticle flux above 7-ray production energies. 

3. (a) The spectrum (as in figure |l]) for explo- 
sion energies C/ = 0.6 x 10^^ , 0.8 x 10^^ , and 
1 X 10^^ erg, all with CO core mass 2.1M0. 
(b) As in (a), with the energies rescaled ac- 
cording to the rule E* = {U/IO^^ eig)~^E; 
for the analytic scaling law this should give 
identical spectra. While this is not per- 
fectly true, clearly the scaling is a good ap- 
proximation. 



4. (a) The spectrum (as in figure ||) for CO 
core masses 1.8, 2.1, and 2.9M0, all hav- 
ing an explosion energy of C/ = 1 x 10^^ 
erg. Although the shapes are similar, the 
higher energy scale for the lower masses 
means that they also have larger intensi- 
ties in the energies above threshold. 

(b) As in (a), with the energies rescaled 
according to the rule E* = {M/2.1Mq)E; 
for the analytic scaling law this should 
give identical spectra and the curves should 
overlap. We see that this scaling is not very 
good, and the results seen more consistent 
with little or no mass dependence. 

5. The yields of LiBeB isotopes func- 
tion of the cutoff energy -Emax- The strong 
-Emax dependence is a result of crossing the 
various thresholds for spallation produc- 
tion. Note that the yields are very small, 
even though our CO core supernova must 
be considered as a best (and uncommon) 
case. Thus it seems unlikely that the LiBeB 
production by fast ejecta is significant for 
chemical evolution. 

6. Elemental and isotopic ratios of LiBeB as 
a function of E^max- The extreme behavior 
at low E^max is due to threshold behavior. 
Since the absolute yields are small, these 
ratios are only of interest if they may be 
observed in a very localized region. The 
strong variations at low E'max are due to 
cross section threshold and resonance ef- 
fects. 
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